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Abstract 

This paper considers the small amplitude unsteady motion of an inviscid non-heat conducting 
compressible fluid on a transversely sheared mean flow. It extends a previous result given in Goldstein 
(1978(b) and 1979(a)) which shows that the hydrodynamic component of the motion is determined by 
two arbitrary convected quantities in the absence of solid surfaces or other external sources. The result is 
important because it can be used to specify appropriate boundary conditions for unsteady surface 
interaction problems on transversely sheared mean flows in the same way that the vortical component of 
the Kovasznay (1953) decomposition is used to specify these conditions for surface interaction problems 
on uniform mean flows. But unlike the Kovasznay (1953) case the arbitrary convected quantities no 
longer bear a simple relation to the physical variables. One purpose of this paper is to derive a formula 
that relates these quantities to the (physically measurable) vorticity and pressure fluctuations in the flow. 

1.0 Introduction 

The small amplitude motion of an inviscid non-heat conducting compressible fluid is governed by the 
Linearized Euler Equations, i.e., the Euler equations linearized about an arbitrary, usually steady, solution 
to those equations — customarily referred to as the base flow. The simplest case occurs when the base 
flow is completely uniform. In a now classical paper, Kovasznay (1953) showed that the unsteady 
isentropic motion on this flow can be decomposed into the sum of a purely vortical disturbance that has 
no pressure fluctuations and an irrotational disturbance that carries the pressure fluctuations. The latter 
satisfies a second-order wave equation when the flow is compressible and should, as argued by Mohring 
(1976), either decay or propagate relative to the base flow. It can, therefore be associated with the 
acoustic component of the motion on these flows. The former, which moves downstream with the mean 
flow, i.e., it is a purely convected quantity, can be associated with the remaining, hydrodynamic, 
component of the motion. Any convected velocity field will satisfy the linearized momentum equation for 
this flow, but continuity only allows two of its components to be arbitrary. These two quantities can then 
be independently specified as steady state boundary conditions for unsteady surface interaction problems. 
This makes the Kovasznay decomposition particularly useful for analyzing the interaction of turbulence 
(which corresponds to the hydrodynamic component of the motion) with surfaces embedded in uniform 
mean flows (Sears, 1941) or in flows that become uniform in the upstream region (Elunt, 1973; Goldstein, 
1978(a), 1979(b)). It is worth noting, however, that the Kovasznay decomposition is not unique because 
there are irrotational (homogeneous) solutions that carry no pressure fluctuations. There have been many 
attempts to extend these ideas to non-uniform base flows, but the situation is considerably more 
complicated when the entire base flow is non-uniform. The simplest case occurs when the base flow shear 
is uniform (i.e., constant). 

In 1907, Orr (1907, pp. 26-29, see Drazin and Reid, 1981, pp. 147-151) pointed out that the linearized 
incompressible vorticity equation 
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where x denotes the time, y,, y 2 - y;; are Cartesian coordinates, with y x being in the mean flow direction or 
equivalently the two-dimensional Rayleigh equation 
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that determine the two-dimensional vorticity perturbation 0)3 and the unsteady transverse velocity 
perturbation v 2 ( vs , x ) for the two-dimensional small amplitude motion on a sheared mean flow with 
linear velocity profile 


U = Xy 2 


(1-3) 


can be integrated to obtain 
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where the transverse vorticity perturbation co 3 is an arbitrary function of its arguments. Orr obtained an 
analytic solution to an initial value problem associated with this equation and used it to study the 
development of the velocity and pressure fluctuations starting from some initial state. But (1.4) can also 
be formulated as a steady-state (i.e., time -stationary) boundary value problem whose formal solution is 
given by 


v 2 ( x ,t) 





CO, 


v Ay 2 j 


dydx 


(1.5) 


where x = {x h x 2 },y= {y \ , >’ 2 } denote the two-dimensional Cartesian coordinates, the inner integration is 
over the entire iq -yq plane, T denotes a large time interval and g 0 denotes a Green’s function, which is 
determined by 


"a 2 | e 2 ' 

8x 2 dx 2 


g 0 (x,t\ y,x) = 5(f-x)5(y-x) 


( 1 - 6 ) 


together with appropriate outgoing wave boundary conditions. 

It can be shown that v 2 (x,t) will decay very slowly with x\ when the arbitrary convected quantity oy 

(x - y 1 fky 2 , y 2 ) decays sufficiently slowly as yq — >■ ±00. The latter can then be specified as an upstream 
boundary condition since(1.5) satisfies (1.2) for any choice of this quantity. Note that the statistics of any 
time -stationary turbulent shear flow can be represented by a sum of such time harmonic functions (see for 
example, Goldstein and Durbin 1980). 

Equation(1.4), which bears some relation to the so-called continuous spectrum, was extended to 
three-dimensional compressible motions on base flows with constant mean shear by Mohring (1976) and 
to general planar compressible shear flows by Goldstein (1978(b), 1979(a)): hereafter referred to as G78 
and G79, respectively) who showed how this more general (non-uniform mean flow) result can used to 
formulate surface interaction problems that are relevant to aircraft noise problems in which the surfaces 
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are typically embedded in transversely sheared mean flows. But, unlike the uniform flow and 
incompressible constant mean shear cases, the arbitrary convected quantity co c (x -y\ /U(y 2 ), yi ) — or its 
equivalent — can no longer be associated with a single component of the vorticity in these more general 
flows. 

G78 and G79 obtained a result that is the natural generalization of the vortical solution in the 
Kovasznay (1953) decomposition by extending the single-convected-quantity, Orr-type, solutions to 
include an additional arbitrary convected quantity, but neither of these can be associated with a single 
component of the vorticity or, for that matter, any other specific physical (measurable) variable. 

Much of modem Rapid Distortion Theory (RDT) is concerned with the small amplitude vortical 
motion on potential flows that originate from a completely uniform upstream flow (Sagaut and Cambon, 
2008; Hunt, 1973, 1977; Goldstein, 1978(a), 1979(b)) or small scale motion on irrotational mean flows 
(e.g., Gartshore, Durbin and Hunt 1983 and Nazerenko, Kevlahan and Drubulle, 1999). An important 
application of existing rapid distortion theory is the specification of appropriate upstream boundary 
conditions, which can only be accomplished by using the Kovasznay decomposition (see, for example, 
Lele 1998, p. 432 and ff. and Colonius 2004, p.316). A natural extension of these ideas would be to 
develop a theory based on the small amplitude vortical motion on steady vortical base flows that originate 
from a parallel, or more generally, a transversely sheared mean flow in the upstream region (Colonius 
2004, pp. 323-324). 

In this paper, we consider the small-amplitude motion of an inviscid non-heat-conducting 
compressible fluid on a general transversely sheared mean flow with the aim of understanding its 
underlying structure and providing the necessary formalism to extend the class of mean flows for which 
applications of RDT can be considered. 

The paper begins by using a newly uncovered connection to the adjoint Rayleigh operator to extend 
the G78 and G79 results to completely general (variable density, non-planar base flow, etc.) transversely 
sheared mean flows. As in G78 and G79 the unsteady motion, which can be calculated by solving an 
inhomogeneous equation, is determined by two arbitrary convected quantities in the absence of solid 
surfaces, external sources and Kelvin-Helmholtz instabilities. The relevant solution is the sum of two 
parts: a part which produces no pressure fluctuations and is driven by one of the two arbitrary convected 
quantities and a part that produces pressure fluctuations and is driven by the other arbitrary convected 
quantity. 

A relation between the convected quantities and the physical variables is worked out in Section 3.0. It 
is shown that there is a specific linear combination of these quantities that is related to a linear 
combination of the physical variables (pressure, vorticity and their derivatives) by a single algebraic 
equation, which generalizes an equation obtained by Mohring (1976) for the linear base flow profile (1.3) 
with constant sound speed. But Mohring (1976) result, which is derived from the Beltrami equation 
(Emanuel, 1993, p. 90), does not include a second convected quantity which, as shown in G79, must be 
included in order to completely represent small amplitude vortical motion on steady vortical base flows. It 
cannot, therefore, be used to represent general three-dimensional vortical disturbances. Our results 
generalize the classical Kovasznay (1953) decomposition to arbitrary transversely sheared mean flows 
and provide the formalism needed to apply RDT to these mean flows. 

The paper is concluded with Section 4.0, which discusses the implication of these results and how 
they can be used to formulate RDT problems, i.e., problems that involve the interaction of turbulence with 
solid surfaces. 

2.0 The Basic Result 

We suppose that the flow is inviscid and non-heat conducting and assume an ideal gas so that the 
entropy is proportional to ln(p/p y ), and the squared sound speed is yp/p, where p denotes the pressure, p 
the density and y the specific heat ratio. Then the inviscid pressure p' = p -p 0 and momentum flux 
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«/ = pvj, 


( 2 . 1 ) 


perturbations (where v' denotes the velocity perturbation)on a transversely sheared mean flow with 
pressure p 0 = constant, velocity U(y T ) and mean sound speed squared c 2 (y T ) are governed by the linearized 
momentum and energy equations 
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is the convective derivative. 

G79 points out that the momentum Equation (2.2) will be identically satisfied for any function <j) and 
any purely convected function 0(x - (y\!U), >’■/) when «, and p' are determined by 
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where 5,, denotes the Kronecker delta, Ujk denotes the alternating tensor and 

A _ d D 0 <\> ! ^ du A> 

J dy j Dx dy j dy l 


(2.7) 


denotes a kind of generalized particle displacement. 

The arbitrary function (j) can then be adjusted to ensure that the energy Equation (2.3) is also satisfied by 
substituting (2.5)— (2.7) into (2.3) to obtain 
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which can then be integrated to obtain 
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where 
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(2.10) 
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It is well known that the momentum flux perturbation, w, can be eliminated between (2.2) and (2.3) to 
show that the pressure fluctuation p' also satisfies Rayleigh’s equation 


Lp' = 0 


( 2 . 11 ) 


where 
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denotes the usual Rayleigh operator. But it is readily verified that 
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for any functions u, v, which shows that the Rayleigh operator/, is adjoint to L a (Morse and Feshbach, 
1953, p. 870). 

Let g(y, x | x,f) 


L g(j,x| x,t) = 8(j-Af)5(x-t) (2.14) 

be a Green’s function for the Rayleigh operator L which exhibits outgoing wave behaviour as | y \ — > go, 
where, as usual, the first two arguments denote the dependent variables and the second two denote the 
source variables. We expect g(y,x \x,t) to vanish as x — > go for all finite y and as vq — > ±co for all finite x 
when the base flow U is globally stable, and causality should ensure that it vanishes as x — > -go. We 
therefore suppose that 


lim g(y,x\x,t), lim g(y,x \ x,t) = 0, (2.15) 

T— »±00 JPj —>±oo 


Setting u equal to g(y,x \ x,t) in (2. 1 3), letting v be a solution to (2.9) and using the divergence theorem 
shows that 
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where T denotes a very large but finite time interval, V is a region of space bounded by cylindrical 
surface(s) S, h = {/?, } is the unit outward-drawn normal to S, where we have omitted terms that are 
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negligibly small as T — > ±co and assumed that the contribution from the end caps can be neglected (Tam 
and Auriault, 1998). This formula expresses the solution to Equation (2.9) in terms of the volume source 
distribution 5) c (x - y\/U(y T ),y T ) and the values of ()> on some arbitrary cylindrical surfaces S. The Green’s 
function g is not uniquely determined by (2.14) and can be required to satisfy certain boundary conditions 
on a portion of the surface S. The present analysis is somewhat unconventional in that the direct Green’s 
function g now plays the role of an adjoint Green’s function for the solution (j). 

When the surface S is at infinity, i.e., V is all space 
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(which generalizes Equation (2.36) in G79) and Equation (2.6) then shows that the transverse velocity 
perturbation v' ± is given by 
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In the more general case, where the surface integral does not vanish. Equation (2.5) shows that the 
pressure perturbation p'(x,t ) is given by 
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when the surface integral vanishes, i.e., when g and (j) satisfy appropriate boundary conditions on any 
cylindrical surfaces that may be present in the flow. 


2,1 Interpretation of Results 

Equation (2. 1 8) can be shown to generalize the Orr result (1.5) by inserting Equation (B. 12) into that 
result, noting that the integral over the second term vanishes and that the Green’s function (1.6) is self- 
adj oint (i.e., go(y,x \ x,t) = go(y,x \ .v,t)) to show that it reduces to ( 1 .5) for two dimensional incompressible 


N ASA/TM— 20 13-217862 


6 



flows with constant mean shear when the arbitrary convected quantity d> c (x -y\/U(y T ), y T ) is replace by 
the renormalized quantity 
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which has dimensions of vorticity (based on the rescaled velocity «,)■ The most significant difference 
between these results is that the convected quantity eo c is no longer equal to the vorticity. 

Mohring (1976) pointed out that the velocity fluctuations on a uniform base flow satisfy a third order 
wave equation which can be integrated to obtain an inhomogeneous second order equation whose 
homogeneous solutions can then be identified with the acoustic motion. But in the present case the 
acoustic motion must be identified with the homogeneous solutions of the third-order wave Equation (2.9) 
, which cannot, in general, be reduced to second order because these solutions have both acoustic and 
hydrodynamic components. It can, however, be reduced to second order when the mean flow does not 
support hydrodynamic motion. This can be done for uniform mean flows by introducing the dependent 
variable ® = D 0 fy/Dx (which can then be used in (2.5)and (2.6) to calculate the pressure and velocity). A 
similar reduction can also be done for two-dimensional flows with constant c 2 and mean shear (see 
Mohring, 1976). This can be seen by differentiating (2.9) with respect to y 2 and using (2.7) and (2.10) to 
obtain 
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which becomes upon using (2.6) to eliminate A 2 
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And using (2.8) to eliminate Dg(|) / Dx 4 leads to Mohring’s (1976) result 
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Mohring (1976) also showed that a similar result can be obtained for the three dimensional motion on 
two-dimensional base flows with constant c 2 and mean shear. 

The unsteady motion is given by (2.5) to (2.7) with 4> determined by the inhomogeneous solution 
(2. 17) of (2.9) for the general transversely sheared mean flow considered in this paper. In the absence of 
solid surfaces and other external sources the resulting unsteady flow consists entirely of subsonically 
propagating disturbances when the flow is subsonic and, therefore, cannot radiate to the far field when the 
flow is unbounded (Goldstein, 2005 and 2009). This can easily be verified in any particular case by 
working out the relevant far field expansion. It is therefore appropriate to identify it with the 
hydrodynamic component of the motion. As in the Kovasznay result this flow is itself the sum of two 
parts, but in this case one of the parts (the part that is driven by the arbitrary convected quantity 
cl),, (x -yJU, >’/)) produces pressure fluctuations while the other part, which is driven by the arbitrary 

convected quantity fl(x —yi/U, >’-/), does not. But unlike the Kovasznay result, the two arbitrary convected 
quantities are not experimentally observable variables. 
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3.0 Relation Between the Physical Variables and the Convected Quantities 


Formulations that represent upstream boundary conditions in terms of the arbitrary convected 
quantities, by (x - yj/U, y T ) and T(x - }j/U, y T ), are most useful when these quantities can be related to the 
physical (preferably measurable) flow variables. This section shows that there is a local (i.e., non-integral) 
equation that relates a linear combination of by (x -y\/U,y T ) and 9(x -y\IU,y T ) to these variables. The 
results are obtained by using relations derived in Appendix A. Relations between the vorticity and 
pressure fluctuations for the combined acoustic-like and discrete hydrodynamic component of the 
unsteady motion are also derived. 

Equations (A.8) and (A. 12) show that 
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which has zero divergence, is the vorticity component generated by the second convected quantity 9, 
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can be interpreted as a transverse particle displacement since it follows from (2.6) that 
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Equations (3.1) and (3.5) relate the arbitrary convected quantities 5> c (x -yi/U(y T ),yr) and 
9(x -yi/U(y T ), yr) to the pressure p', the transverse vorticity components s ;//c (OjN k and the particle 
displacement while Equations (3.2) and (3.5) relate the arbitrary convected quantity 9(x -y\/U(y T ), yr) 
to the pressure p',the normal vorticity component o OjNj and the particle displacement q i . But the particle 
displacement q ± (which, unlike the pressure and vorticity components, is not a physical variable in the 
usual sense) can be eliminated between Equations (3.1) and (3.2) when dNJdyj ± 0 to obtain the following 
equation 
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of the two arbitrary convected quantities cb c (x - y\/U(y T ),y T ) and 9(x -y\/U(y T ),y T ) to the physical 

variables p’ and co,. Equations (3.3), (3.10), (3.1 1) and (A. 7) show that % c (x -y\/U(y T ), yr) is itself an 
arbitrary convected quantity. 


3,1 Application to Planar Mean Flows 

These results are fairly complex, but most of the relevant literature is concerned with the small 
amplitude motion on planar base flows, where c 2 and U depend on a single Cartesian coordinate (sayy 2 ). 
Equations (3.1), (3.2), and (3.5) then become much simpler and reduce to 
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(3.14) 


where 


_d_ 

8y\ 


a> c -p — 


dN ? 
dy : 


-hi 


= N, 


dco 3 

dy x 


Scoj ' 

A, 


— V,?d = <B, + 
N 2 11 



V 


2 

II 


a 2 | d 2 

ty\ dy\ 


(3.15) 


(3.16) 


These equations show that the arbitrary convected quantities 6) c (x - y\/U(y T ), y T ) and T(x -y\/U(y T ), y T ) 

are now related to p' the vorticity components coi, co 3 and the particle displacement q by fairly simple 
relations. Equations (3.12) and (3.13) now simplify to 


Xc(^-y^u(y 2 ),y T ) 


d dp' 
dy x \ L dy 3 





where 


lc(i-y\IU{y 2 ),y T ) 


d f d(b c 

l A 


N 2 dy 2 


(3.17) 


(3.18) 


and we have used the fact that the solenoidal vectors co, and co-* 1 have zero divergence. 

Notice that Equation (3.14) is now independent of ooj°° 1 and, therefore of the second convected 

quantity 9. But the divergence dNJdy, is equal to zero for the constant shear-constant c 2 parallel mean 
flow (1.3), since TV) is a is constant in that case and it follows from (2.22) that Equation (3.14) then 
reduces to Mohring’s (1976) result 


A" 

fyl V 


cm 



c j 


dco 3 oc»)| 
dy x dy 3 


(3.19) 


that relates the single convected variable eo c to the physical variables. Equation (3.17) relates the 
combined convected quantity ~/ c ( T ~ydU(y 2 ), Jr) to the physical variables // and co,. 


4.0 Concluding Remarks 

This paper is based on a formal solution to the linearized Euler equations for a transversely sheared 
mean flow which shows that the general solution to these equations is given by an inhomogeneous 
solution to the adjoint Equation (2.9) and, like the Kovasznay decomposition for the unsteady motion on 
uniform flows, involves two arbitrary convected quantities T(x -y x IU,y T ) and cb f (x - y\IU,y T ) that can 
be used to specify appropriate boundary conditions for unsteady surface interaction problems and, 
thereby, extend the class of flows for which this can be done to arbitrary transversely sheared mean flows. 
It is shown that this equation reduces to the second order wave equation derived by Mohring (1976) for 
two dimensional motions on the linear base flow (1.3) with constant sound speed once an appropriate 
dependent variable is introduced. But unlike the Kovasznay result, the hydrodynamic component carries 
pressure fluctuations and the acoustic component is not irrotational. Section 3.0 shows that there is a 
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linear combination of the convected quantities that is equal to a linear combination of the pressure and 
vorticity 

The results developed in this paper are expected to be particularly useful for representing incident 
vortical disturbances in Rapid Distortion Theory problems involving the interaction of turbulence with 
surfaces embedded in transversely sheared base flows or, more generally, in vortical base flows that 
asymptote to transversely sheared mean flows in the upstream region. Equations (2.18), (2.21), and 
similar equations that can be derived for the remaining velocity components, determine the velocity and 
pressure fluctuations in terms of the arbitrary convected quantities by (x - y\IU,y T ) and 0(x - y\/U,y T ) in 
these flows. 

One such application currently under investigation is to the noise generated by the interaction of a 
high Reynolds number rectangular jet with a nearby flat plate (Afsar, Goldstein and Leib 2013). This 
work applies the formalism developed in the present paper to predict the acoustic spectrum of the noise 
scattered from the trailing edge of the plate and compares it with experimental data. The predictions 
require a model for the autocovariance of the convected quantity by (x -y\/U, y T ), whose validity can be 
tested by using the relations developed in Section 3.0 to relate it to physical quantities that can be 
measured in the laboratory. 
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Appendix A. — Relation Between Convected Quantities and Physical Variables 


Equations (2.5) to (2.7), (2.9)and (2.10) imply that the convected quantity eb c is related to the pressure 
perturbation p' and the generalized particle displacement A, by 


, 8 2 
=P +—c A i 
8y t 


D, 


u i “TT" A ; _ S n r l_ L +£ 


1 8U 8 


Dx ' 11 ijk c 2 dy, 8y k V t/ 


■a 


x_ 77 , d ; r 


(A.l) 

(A.2) 


where is r| ± given by (3.6). 

The density-weighted vorticity eo, is therefore, given by 

8u k 8A k 8U 8A k Sri, 

= S 'A "^7” “ 7>r S,/Vt ■ + S ' Jk ^ ^ S ' 71 + &knm& v k 

But since Equation (2.7) and the identities 


8yj 8y x 8y j 


8 i su aa 
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show that 


this can be written as 


d 2 U 


J yk 8y jdy k &ljk 8yjdy k 
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= 0 
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Ijk 8y j yk 8y k 8y x 8y 
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(A. 3) 


(A.4) 


(A. 5) 
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Dx 8y t 


+ co) 


k J 


where co-*' is given by (3.5). The definition (2.7) and the identity 

8U 8U 


e ijk ‘ 


tyj fyk 

can now be used to show that (A.6) can be written as 


- = 0 


CO, = -8 


3r|j 


— p 

ii\ ^ ijk ^ 

dvj oy j 8y x 
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k +®!” ) 


(A.6) 


(A. 7) 


(A.8) 


And it, therefore, follows from (3.6) that the cross product s^eo jdU / 8y k can be written as 
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which can be solved for 5A ; / 5a) to obtain 
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Inserting this into (A. 1 ) shows that 
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where the identity (A. 7) has again been used. 


(A. 9) 


(A. 10) 


(A- 11) 


(A. 12) 


(A- 13) 
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Appendix B. — Limiting Form of Equation (2.18) for Linear Base Flows 


The easiest way to show that Equation (2. 1 8) reduces to Equation (1 .5) is by noting that the Fourier 
transform 


_ 1 r r /,x/ x / 
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which satisfies the reduced Rayleigh equation 
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since g(y, x | x,t) depends on x, t and y u x\ only in the combinations t-x and x\ -y\, respectively. 
It, therefore, follows that 
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for planar base flows with constant sound speed — which means that 
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(B.7) 
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for the linear profile (1.3). And it follows from (B.8) that 
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Taking the inverse Fourier transform of (B.9) shows that 


g 2 ( y,T | x,t) = {y, t I x,t) + ( JA I x,t ), 


- 2 dx x 
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while the limiting form of Equation(B.l 1) for a completely two dimension flow is just the Fourier 
transform of Equation (1.6). 
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